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ABSTRACT 


i  i  .  \  *_-■  <  i  >  .  ■■ '  I 

The  dissipation  method  is  used  to  obtain  estimates  for 
the  friction  velocity,  U*,  as  well  as  values  for  the  neutral 
drag  coefficient,  C^N,  for  data  collected  from  a  coastal 
tower  off  San  Diego,  California.  CDI^  is  found  to  be  indepen¬ 
dent  of  the  ten-meter  height  windspeed,  U-^q,,  for  velocities 
between  4-9  m/sec.  Its  value  is  estimated  to  be 
(0 . 94+0 . 4) 10J  which  compares  well  with  values  by  Smith 
(1980)  and  Large  and  Pond  (1981).  Definite  trends  in  CD^ 
with  fetch  and  sea  state  are  also  observed.  Drag  coeffcient 
estimates  are  found  to  be  higher  for  short  fetch  than  for 
long  fetch  conditions.  is  also  seen  to  increase  sharply 

just  before  frontal  passages  and  during  sea  breeze  condi¬ 
tions  when  the  waves  are  actively  growing.  With  the  wind- 
speed  and  wave  field  reaching  equilibrium,  CDN  is  found  to 
decrease  with  time  to  a  smaller  and  more  constant  value. 
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I.  INTRODUCTION 

The  Marine  Atmospheric  Boundary  Layer  (MABL)  is  the 
region  adjacent  to  the  air- sea  interface  where  the  vertical 
variation  of  the  mean  stresses  is  large  compared  to  the 
horizontal  gradients  (Holton,  1979).  In  this  layer,  turbu¬ 
lent  shear  stresses  due  to  both  molecular  properties 
(viscosity  and  diffusion)  and  turbulent  eddies  coexist. 
These  eddies  are  very  effective  mixing  agents  which  serve  to 
transfer  momentum,  heat  and  other  quantities  to  and  from  the 
earth’s  surface,  at  a  rate  much  faster  than  that  for  molec¬ 
ular  diffusion.  Frictional  forces  due  to  molecular  viscosity 
are  usually  neglected  by  means  of  scale  analysis  (Holton, 
1979),  and  the  emphasis  has  been  placed  on  parameterizing 
the  fluxes  resulting  from  the  turbulent  eddies. 

The  mean  vertical  fluxes  of  momentum,  sensible  heat  and 
latent  heat  are  (Fleagle  and  Businger,  1980): 


"Pau  w 
Cpfl'w' 

Lfw ' q ' 


Ua) 

(lb) 

(lc) 


where  r,  H  and  E  are  the  momentum  flux  (wind  stress), 
sensible  heat  flux  and  latent  heat  flux,  respectively.  u' 
and  w'  are  the  horizontal  and  vertical  turbulent  wind  compo¬ 
nents,  while  Q'  and  q'  represent  the  turbulent  fluctuations 
of  potential  temperature  and  specific  humidity.  Lf  is  the 
latent  heat  of  vaporization,  pa  is  the  air  density  and  Cp  is 
the  specific  heat  of  air  at  constant  pressure. 

In  order  to  estimate  these  fluxes  using  routine  surface 
layer  observations,  bulk  transfer  formulae  were  developed 
(see  Fleagle  and  Businger,  1980): 


PaCD<Uz  •  Uo)' 


(2a) 
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Pa^pCH'Uz  "  Uo)(To  -  V 
PaLfCE(Uz  -  U0)(qQ  -  qz) 


(2 
(2c) 


where  U2 ,  Tz  and  qz  represent  the  windspeed,  temperature  and 
humidity  specified,  by  convention,  for  a  reference  height, 
z.  UQ,  TQ  and  qQ  are  the  same  quantities  measured  at  the 
surface.  Cq,  and  Cg  are  dimensionless  quantities  repre¬ 
senting  bulk  transfer  coefficients  for  momentum,  sensible 
heat  and  latent  heat,  and  are  for  the  height,  z.  These  are 
known  as  the  drag  coefficient,  the  Stanton  Number,  and  the 
Dalton  Number,  respectively. 

The  turbulent  exchange  processes  greatly  influence  the 
general  circulation  of  the  atmosphere  and  thus  have  tremen¬ 
dous  impact  on  countless  social  and  environmental  concerns. 
Central  to  their  understanding  is  a  determination  of  the 
nature  of  the  bulk  transfer  coefficients,  particularly  the 
drag  coefficient,  Cq.  This  quantity  encompasses  the  mecha¬ 
nism  responsible  for  the  wind  stress  at  the  surface.  Since 
most  of  the  parameters  used  to  describe  the  bulk  formulae 
are  routinely  measured,  ascertaining  a  reliable  value  for  CD 
based  on  these  quantities  is  essential  for  accurately 
defining  the  surface  stress,  t.  Further,  CD  is  an  important 
factor  driving  many  of  the  current  ocean  wave,  acoustic  and 
boundary  layer  models  (Geernaert,  1985a). 

Considerable  effort  has  gone  into  examining  the  over¬ 
water  drag  coefficient.  It  has  been  shown  to  be  windspeed 
and  stability  dependent,  and  may  also  vary  with  factors  such 
as  fetch  and  depth  (Geernaert,  1983).  Studies  also  indicate 
that  the  overwater  drag  coefficient  tends  to  have  a  magni¬ 
tude  on  the  order  of  10'^.  Its  value  exhibits  considerable 
scatter  under  varying  environmental  conditions.  Higher 
values  are  usually  associated  with  an  unstable  atmosphere, 
short  fetch  and/or  shallow  seas,  implying  a  wave  dependence 
(Geernaert,  1983).  The  more  recent  use  of  remote  sensing  to 


12 


study  winds  over  the  ocean  requires  an  understanding  of  this 
wave  dependence  since  it  assumes  that  the  wave  structure  is 
related  to  the  wind  stress.  Clearly,  the  processes 
affecting  the  drag  coefficient  merit  considerable  further 
study.  This  thesis  will  investigate  the  dependence  of  the 
drag  coefficient  on  windspeed,  fetch  and  wave  state. 

An  important  scaling  parameter,  based  on  the  wind 
stress,  is  the  friction  velocity  ,  U*,  which  may  be  defined 
as : 

T  =  paU*2  (3) 

The  quantity  U*  is  related  to  most  air- sea  exchange 
processes  such  as  heat  and  moisture  fluxes  (Businger,  1973), 
the  surface  drift  (Hicks,  1972)  and  ocean  mixed  layer 
dynamics  (Kraus,  1972). 

Data  considered  in  this  investigation  were  from  the 
Tower  Ocean  Wave  and  Radar  Dependence  (TOWARD)  experiment. 
Located  offshore  of  Mission  Beach  at  San  Diego,  California, 
the  Naval  Ocean  Systems  Center  (NOSC)  tower  provided  a 
stable  platform  with  a  capability  for  mounting  a  variety  of 
sensors  to  measure  the  variables  needed  to  parameterize  the 
bulk  formulae.  Naval  Postgraduate  School  (NPS)  personnel 
made  atmospheric  surface  layer  measurements  of  both  mean  and 
fluctuating  velocity,  temperature  and  humidity.  The  NPS 
measurements  are  listed  in  Table  I.  Wave  information  based 
on  radar  signal  returns  was  collected  through  investigations 
from  the  Naval  Research  Laboratory,  the  Jet  Propulsion 
Laboratory  and  the  University  of  Kansas.  The  Scripps 
Institute  of  Oceanography  measured  tower  frequency  wave 
spectra  with  an  array  of  sub-surface  pressure  transducers. 

The  measurement  strategy  for  the  TOWARD  Experiment  was 
originally  composed  of  two  intensive  data  collection 
periods:  September/October  1984  and  January /  February  1985. 
These  periods  were  chosen  because  of  their  differing  average 


windspeeds  and  resultant  sea  states.  With  its  relatively  low 
average  windspeeds,  the  fall  measurement  period  had  calm  to 
moderate  sea  state  conditions.  The  winter  period,  with  its 
higher  average  windspeeds,  was  designed  to  focus  on  more 
active  sea  states.  Equipment  problems,  however,  forced  the 
postponement  of  Phase  I  to  October/November  1984  and  Phase 
II  to  March/April  1985. 

Figs.  1.1  and  1.2  provide  two  views  of  the  NOSC  tower. 
Located  in  17-18  meters  water  depth,  the  tower  provided  a 
stable  platform  for  data  collection.  A  structural  extension 
of  the  tower  from  the  south  side  allowed  wave  rider  sensors 
to  be  placed  away  from  the  main  tower  structure. 

Measurements  of  wind  speed  and  direction,  wave  height, 
tide  level  and  near  surface  currents  were  recorded  routinely 
during  intensive  measurement  periods.  Air  and  sea- surface 
temperatures  were  also  recorded,  along  with  wind  stress  data 
that  was  obtained  with  a  sonic  anemometer  and  a  hot  film 
sensor . 
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TABLE  I 

NPS  Measurements 


Measurement 

Sea- Surface  Temperature 

Mean  Surface  Layer: 

Wind  (speed  and  direction) 

Temperature 

Humidity 

Radiation  (solar  and  IR) 

Turbulent  Surface  Layer: 

Wind  stress:  direct 
Wind  stress:  indirect 
K.  E.  dissipation  rates 
Sensible  heat  flux 

Latent  heat  flux 

Spectra:  windspeed 
temperature 
humidity 


Sensor 

Floating  platinum 
thermometer 


Cup  anemometer, 
vane 

Aspirated  platinum 
thermometer 
Dew  cell 
(cooled  mirror) 
Eppleys  (up  and 
down) ,  Total 


Sonic  anemometer 
Hot  films 
Hot  films 

Platinum  thermometer 
and  sonic  anemometer 
Lyman- alpha  and 
sonic  anemometer 
sonic  anemometer 
Platinum  thermometer 
Lyman- alpha 


Figure  1.2  The  NOSC  Tower 


II.  REVIEW  OF  EXISTING  SURFACE  LAYER  THEORY 


The  surface  layer  comprises  approximately  the  lowest  ten 
per  cent  of  the  MABL  and  is  defined  as  the  layer  of  air 
nearest  the  ground  where  atmospheric  parameters  may  be 
scaled  with  surface  fluxes.  Within  this  layer,  an  approxima¬ 
tion  adopted  is  that  the  fluxes  are  independent  of  height 
(Businger,  1973). 

In  Eqn.  (2a),  the  surface  speed,  UQ ,  is  a  parameter 
necessary  for  estimation  of  the  drag  coefficient.  However, 
because  UQ  is  very  small  compared  to  U-^q,  Uq  is  assumed  to 
be  zero  in  the  bulk  definitions.  With  this  assumption,  Eqns . 
(2a)  and  (3)  may  be  combined  into  the  form: 

CD10  =  (U*/U1Q)2  (4) 

Note  that  the  drag  coefficient  in  Eqn.  (4)  is  defined  for  a 
height  of  10  meters  and  is  denoted  as  Cp^g. 

The  similarity  theory  which  described  the  mean  surface 
layer  wind  profile  is  (Monin  and  Obukhov,  1954): 

dU/dz  =  (U*/kz)<£m(z/L)  (5) 

where  dM/Qz  represents  the  vertical  gradient  of  the  wind 
speed  and  k  is  the  von  Kartnan  constant,  which  ranges  in 

value  from  .35  to  .43  (Panofsky  and  Dutton,  1984).  This 
study  assumes  a  value  of  0.4  for  this  constant.  <£m(z/L)  is 
a  dimensionless  stability  function.  One  parameterization  of 
this  function  is  that  presented  by  Businger,  et .  al . , 
(1971) : 

<>m  =  [1  -  a(z/L)]-1/4  ( z/L<0 )  (6a) 

<f>m  =  [1  +  J8(z/L)]  ( z/L>0 )  (6b) 
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where  a  and  j3  are  diabatic  profile  constants,  with  approxi¬ 
mate  values  of  15±1  and  5i0.5,  respectively  (Geernaert, 
1983).  z  is  the  measurement  height  and  L  is  the 
Monin-Obukhov  length,  defined  as  the  height  at  which  the 
buoyant  energy  production  equals  the  shear  production 
(Businger,  1973).  z/L  is  defined  as: 

z/L  =  -gkzw'V/V  (U*3)  (7) 

The  virtual  temperature  fluctuation,  Tv' ,  is  related  to 
the  air  temperature  fluctuation,  T',  and  specific  humidity 
fluctuation,  q',  by: 

Tv'  =  T' (1  +  . 61q)  +  . 61Tq'  (8) 
Thus : 

W7 =  vTT’  (1  ♦  .  6  lq )  +  .eiTvTq'  (9) 

where  g  is  the  gravitational  acceleration.  The  first  term 
on  the  right  hand  side  of  Eqn.  (9)  represents  the  tempera¬ 
ture  flux  and  the  second  term  describes  the  humidity  flux. 
These  terms  are  further  estimated  by  bulk  equations  where: 


“•1"  =  ChU10(To  -  Tz)  (10a) 

w'q'  =  CEU10(qo  -  q2 )  (10b) 

Both  Cg  and  Cj|  are  on  the  order  of  10"^,  but  are  less  well 
defined  than  the  drag  coefficient  (Geernaert,  1983). 

Under  neutral  or  adiabatic  conditions, 

<£m(z/L)  =  1  (11) 

and  Eqn.  (5)  reduces  to: 

dU/dz  =  U*/kz  (12) 

Integrating  Eqn.  (12): 

■CJ>  dU  =  (U*/k)/J^  dz/ 2  (13) 


yields  the  well  known  logarithmic  wind  profile: 


u 


(U*/k)ln(z/z0) 


(14) 


where  the  roughness  length,  zQ ,  is  a  measure  of  the  effec¬ 
tive  surface  roughness. 

For  non-neutral  (diabatic)  conditions,  the  stability 
parameter  must  be  included.  Integrating  Eqn.  (5)  yields  the 
stability  dependent  logarithmic  wind  profile  (Paulson, 
1970) : 

U  =  (U*/k) [ln(z/zQ)  -  xp]  (15) 

where  xp  represents  the  stability  parameter. 

For  unstable  stratifications  (z/L<0): 

xp  =  21n[ ( 1+x) / 2]  ♦  ln[(l+x2)/2]  -  2(tan)‘1x  +  (n/2)  (16) 

where : 

x  =  [(1  -  aU/L)]1/4  (16a) 

For  stable  conditions  (z/L>0): 

xp  =  -  0(z/L)  (17) 

The  logarithmic  wind  profile  is  important  primarily  for 
two  reasons.  First,  it  can  be  used  to  estimate  the  ten 
meter  height  wind  speed,  U-^q.  Beginning  with  Eqn.  (14),  the 
wind  speed  at  a  particular  height,  z,  may  be  described  as: 

Uz  =  (U*/k)  [ln(z/z0)  -  xpz]  (18) 

At  ten  meters  elevation: 

U10  =  (U*)/k[ln(10/zo)  -  ^10]  (19) 

Combining  Eqns .  (18)  and  (19)  yields  the  following  relation¬ 
ship  : 

U10  =  U2  *  (U*/k)[ln(10/z)  -  xpz  ♦  ^10]  (20) 

This  particular  equation  is  important  because  it  is  not 
always  possible  to  place  sensors  on  ships  or  platforms  at 


the  required  ten  meter  height.  As  a  result,  U^q  can  be 
estimated  regardless  of  the  height  of  the  instrument. 

Second,  the  log  wind  profile  is  important  in  determining 
the  neutral  drag  coefficient,  C^jj,  which  may  be  defined  as 
the  wind  stress  exerted  on  the  surface  for  neutral  stratifi¬ 
cations.  Estimation  of  the  neutral  drag  coefficient  is 
important  because  it  eliminates  variations  in  Cp  parameteri- 
zations  due  to  differing  stability  conditions.  Many 
modellers  prefer  to  be  provided  with  C^N  instead  of  CD 
because  the  in-situ  bulk  weather  data  may  be  used  to  esti¬ 
mate  a  stability  correction.  Eqn.  (24)  may  then  be  applied 
to  a  given  CjjN  and  i p. 

A.  GENERAL  DRAG  COEFFICIENT  REPRESENTATIONS 

Eqns .  (4)  and  (15)  yield: 

CD  =  (U*/U2)2[k/[ln(z/z0)  -  ^m(z/L)]]2  (21) 

Under  neutral  conditions,  where  \p  =  0: 

CDN  =  [k/ln(z/zon)J2  (22) 

where  CDN  and  zQn  are  the  drag  coefficient  and  roughness 
length  under  neutral  conditions.  Rearranging  Eqn.  (22) 
yields : 

kcDN'1/2  =  lnU/zon>  (23) 

Combining  Eqns.  (21)  and  (23)  produces  the  following  general 
result  (Geernaert,  1983): 

CD  =  [CDN'1/2  -  \p/k  -  (l/k)ln(z0n/z0)]'2  (24) 

The  drag  coefficient  is  known  to  vary  under  differing 
environmental  conditions.  Studies  have  indicated  that  CD 
increases  more  rapidly  with  increasing  winds  over  water  than 
over  land  (Hsu,  1974).  Geerneart  (1983)  discussed  the 

application  of  Eqn.  (24)  to  both  land  and  water  cases.  Over 


land,  the  roughness  length  is  assumed  to  be  largely  indepen¬ 
dent  of  the  surface  stress.  Thus,  zon  =  zQ  and  Eqn.  (24)  is 
reduced  to: 

CD  =  (CDN_1/2  -  */k)"2  (25) 

Over  water,  both  the  drag  coefficient  and,  hence,  the 
roughness  length  show  considerable  windspeed  and  stress 
dependence  (Kitaigorodskii ,  1973).  Chamock  (1955)  was 
among  the  first  to  recognize  that  increases  considerably 
with  windspeed  over  the  sea.  He  postulated  that  the  rough¬ 
ness  length,  and  hence,  the  neutral  drag  coefficient  was 
proportional  to  the  wind  stress  through  the  following 
scaling  argument: 

z0  =  a'U*2/g  (26) 

where  a’  is  the  Charnock  constant  with  a  typical  value  of 
.015-. 020  (Geernaert,  1983).  Charnock' s  work  served  as  a 
basis  for  subsequent  parameterizations  of  the  drag  coeffi¬ 
cient  and  the  roughness  length  in  terms  of  sea  state 
elements . 

Kitaigorodskii  (1973)  further  described  zQ  based  on  wave 
spectral  densities.  Since  that  time,  zQ  and  CD  have  been 
shown  to  be  dependent  on  a  variety  of  other  factors  that  are 
linked  to  the  character  of  surface  waves  such  as  fetch, 
depth  and  wave  age.  Cp,  therefore,  can  be  evaluated  as: 

CD  =  CDN(U1Q,  z/L,  waves)  (27) 

As  a  result,  a  different  relationship  is  often  used  to  esti¬ 
mate  the  drag  coefficient  over  the  ocean.  Applying 
Charnock’ s  relation  to  Eqn.  (24)  yields  (Geernaert,  1983): 

CD  =  [CDN'1/2  -  */k  *  (l/k)ln(CDN/CD)]'2  (28) 
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B.  DETERMINING  THE  WIND  STRESS  BY  THE  DISSIPATION  METHOD 


Since  estimation  of  the  drag  coefficient  requires  values 
for  the  friction  velocity,  it  is  important  to  examine  how 
the  momentum  flux  is  measured.  Two  primary  methods  for 
obtaining  j  are  the  eddy  correlation  technique  and  the 
dissipation  method.  The  former  method  involves  measurement 
of  the  turbulent  wind  components  u' ,  v'  and  w' ,  and  ideally 
requires  a  steady  platform.  Because  of  the  low  frequency 
motions  inherent  in  sea-going  vessels  on  which  most  marine 
studies  are  conducted,  the  dissipation  method  is  often 
substituted  for  oceanic  wind  stress  experiments. 

The  dissipation  method  is  an  indirect  approach  in  which 
the  rate  of  dissipation  of  turbulent  kinetic  energy  is  used 
to  determine  U*.  This  method  requires  application  of  the 
turbulent  kinetic  equation: 

de/dt  =  U*2dU/dz+ (qw'Tv' )/T0-d/dz[w' e+ (l/pa)w'p ' ] -e  (29) 

where  de/dt  is  the  local  time  rate  of  change  of  turbulent 
kinetic  energy  (TKE),  U*2dU/dz  is  the  wind  shear  produc¬ 
tion,  (qw'Tv«)/T0  represents  the  buoyancy  forces, 
d/dz[w'e  +  (l/pa)w'p']  is  the  redistribution  of  the  flux  of 
TKE  and  work  done  by  pressure  fluctuations  and  £  is  epsilon, 
the  rate  of  dissipation  of  turbulent  kinetic  energy. 

Steady  state  conditions  are  assumed.  The  divergence  term 
is  small  relative  to  the  remaining  terms  in  the  equation  and 
is  therefore  neglected  (Large  and  Pond,  1981).  Eqn.  (29) 
then  simplifies  to: 

U*2dU/dz  +  C(g^V)/T0]  -  £  =  0  (30) 

After  combining  Eqns .  (5),  (7)  and  (30),  and  rear¬ 
ranging,  an  expression  for  the  friction  velocity  is 
obtained : 


1/3 


[ekz/[<£  (z/L)  -  (z/L)]] 


(31) 


For  unstable  conditions,  the  denominator  is  often  replaced 
by: 

<f>£( z/L)  =  0m(z/L)  -  (z/L)  (32) 

Eqn.  (31)  may  then  be  written  as: 

U*  =  [£kz/<^e(z/L)]1/3  (33) 

with  stability  corrections  for  dissipation  scaling  (Wyngaard 
and  Cote ' ,  1971) : 

<f)e  =  (1  +  0.51  z/L|  3/5)3/2  z/L<0  (34a) 

<fie  =  (1  +  2 . 5 1  z/L( 2/3  )3/2  z/L>0  (34b) 

The  friction  velocity,  U*,  can  therefore  be  estimated  in 
terms  of  the  dissipation  rate,  £.  Kolmogoroff  theory  is  then 
employed  to  determine  a  value  for  £  based  on  high  frequency 
wind  fluctuations.  All  spectra  can  be  divided  into  three 
portions:  the  energy  producing  subrange,  the  inertial 

subrange  and  the  dissipation  subrange  (Panofsky  and  Dutton, 
1984).  Large  eddies  from  the  energy  containing  subrange 
cascade  kinetic  energy  and  momentum  down  to  smaller  eddies, 
eventually  reaching  the  dissipation  region  where  molecular 
viscosity  converts  the  kinetic  energy  to  heat.  Between  these 
two  regions  is  the  inertial  subrange.  In  this  region,  energy 
is  neither  produced  nor  dissipated.  Kolmogoroff  theory 
states  that  in  the  inertial  subrange  of  isotropic  turbu¬ 

lence,  the  high  frequency  fluctuations  of  the  wind  velocity, 
u'  can  be  expressed  by  a  one-dimensional  power  spectrum 
(Schacher,  et .  al.,  1981): 

Su(k*)  =  a*€2/3k*'5/3  (35) 

where  Su  is  the  energy  density  spectrum  of  the  horizontal 

windspeed,  is  a  coefficient  of  magnitude  0.50  (Champagne, 

et.  al.,  1977),  and  k*  is  the  wavenumber. 


Applying  Taylor's  frozen  turbulence  hypothesis: 


k*  =  2irf/U  (36) 

equation  (35)  becomes: 

Su(k*)  =  a*€2/3(2fl-f/U)"5/3  (37) 


The  Kolmogoroff  Spectrum  can  now  be  integrated  between 
two  wavenumbers,  designated  k*i  and  k*^  to  give  the  energy, 
a^Ak-'-  interval  (Schacher,  et.  al.,  1981): 
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a  Ak*  = 

j/T/  Su(k*)dk* 

(38) 

k. 

:>h'2/3) 

(38a) 

Combining  Eqns . 

(36)  and  (38)  yields 

: 

2 

°  Af 

0.230(£U)2/3(f1'2/3 

- 

v2/3> 

(39) 

where  f1=(k*1U) 

/ 2  and  fh=(k*hU)/2. 

Rearranging  Eqn. 

(39)  by 

using  Eqn .  (35) 

produces : 

U*  =  2.6  5aU(^fCz 

/ [U0e(z/L)]] 1‘  3( f l~2 

3 

-  -V2  V  2 

(40) 

For  the  TOWARD 

Experiment,  =  5 

Hz 

and  =  50 

Hz ,  and 

z=22  meters.  Su 

bstitutmg  these  valu 

es 

reduces  Eqn. 

( 40  )  to : 

U* 

7  ■  81<7UAf  17  L  '  ■* 

1 

3 

(41) 

Because  <££(z/L 

)  also  depends  on 

Eqn .  (41) 

must  be 

solved  iteratively. 


C.  IN- SITU  HOT  FILM  CALIBRATION 


Calibration  of  the  hot  film  sensor  was  accomplished  by 
recording  the  changes  in  voltage  produced  by  the  winds 
passing  over  the  sensor.  The  relationship  may  be  expressed 
by: 


V2 


1/2 


♦  BUrel 


(42) 


where  V  is  the  voltage,  VQ  is  the  plotted  calibration  inter¬ 
cept  and  tTrel  is  the  relative  windspeed.  B  represents  a 
calibration  factor  that  must  be  continuously  determined  for 
each  data  collection  period. 

B  may  be  estimated  either  through  the  dynamic  or  the 
regression  methods.  The  dynamic  method  begins  by  differen¬ 
tiating  Eqn.  (42): 


4VUrel1/2(dV/dU) 


(43) 


Over  each  data  collection  period,  dV/dU  is  estimated  by 
av/aU’  where  ct\j  and  ay  are  the  standard  deviations  of  U  and 
V,  respectively.  Averaging  V  and  Orel  and  solving  for  B 
yields  : 


4VUrel1/  2  (ay/  a{j) 


(44) 


Rearranging  Eqn.  (44)  and  substituting  into  Eqn.  (40) 
produces  the  following  estimate  for  U* : 


2.81[(4V/B)Urel1/2]av4fCz/U^e(2/L)]1/3(f1'2/3  -  fh'2/3)  (45) 

The  regression  method  is  based  on  values  of  V  and  Urel 
over  several  time  periods  with  varying  average  relative 
windspeeds .  V  is  then  usually  plotted  as  a  function  of  Urel. 
According  to  Eqn.  (42),  the  points  should  fall  along  a  line 
with  a  slope,  B,  where  VQ  represents  the  y-intercept. 

Measurement  -  relative  scatter  in  both  Urel  and  V  can 
cause  errors  in  the  slope.  Therefore,  it  is  sometimes  advan¬ 
tageous  to  choose  V  as  the  independent  variable  if  it  shows 
less  scatter.  In  cases  where  it  is  difficult  to  determine 
which  variable  is  the  least  reliable,  two  regressions  are 
calculated,  with  Urel  and  V  alternated  as  independent  vari¬ 
ables.  B  is  then  determined  from  the  average  of  the  two 


D.  WIND- WAVE  COUPLING  DYNAMICS 

The  total  wave  field  energy  is  divided  into  swell  and 
wind  wave  components: 


=  a 


2  2 
_  +  IT  _ 


(46) 


where  represents  the  total  energy  in  the  wave  field 

2  2 

and  or  g  and  o  ^  are  the  energies  contained  by  the  swell  and 

O 

wind  waves,  respectively,  a  ^  cannot  be  predicted  based  on 
local  scaling,  a  w  has  received  considerable  study  and  a 
great  deal  of  theory  exists  concerning  this  parameter.  The 
approach  currently  is  to  describe  a  in  terms  of  a  local 
equilibrium  spectrum.  This  theory  assumes  that  wind  forcing 
controls  the  dynamics  of  every  wave  component  of  the  wind 
wave  spectrum.  Each  wave  is  characterized  by  an  angular 
frequency,  n' ,  where: 

n’  =  2nf  (47) 

Gravity  is  the  primary  restoring  force  for  these  waves. 

The  local  equilibrium  energy  density  spectrum,  S(n'),  is 
defined  as: 


/T  S  (n '  )dn' 


(48) 


and  S(n')  may  be  parameterized  in  terms  of  Phillips  scaling 
arguments  (Phillips,  1980): 


S(n’) 


j3g(n')' 


n  >n. 


(49) 


Beta  is  the  Phillips  coefficient  which  can  vary  from  .006  to 
.020  (Geernaert,  1983).  nQ '  represents  both  the  wind-wave 
peak  frequency  and  the  boundary  between  the  swell  and  equi¬ 
librium  ranges.  Waves  of  frequency  n'<n0'  represent  swell. 
Wind  waves  are  defined  to  be  those  with  frequencies  of 
n’>n0’.  nQ '  is  usually  approximated  by: 


g/UiQ 


(50) 


for  steady  state  conditions  (Kraus,  1972). 

Beta  is  often  described  in  terms  of  dimensionless  fetch, 


gx/(U*2) 


(51) 


where  x  is  the  upwind  fetch.  Large  values  of  are  associ¬ 
ated  with  small  values  of  fetch.  Beta  has  also  been  parame¬ 
terized  over  an  applicable  range  of  .25  to  1.0  hz  by 
(Geemaert,  et.al.,  1985a): 


0.005  ♦  . 002n'  ♦  1.5(Cn/U*)"2 


(52) 


where  CQ  represents  the  phase  speed  of  the  largest  wind 
wave.  (CQ/U*)  is  defined  as  the  wave  saturation  parameter, 
and  this  parameter  behaves  in  the  same  manner  as  the  wave 
age  C/U.  (C  /U*)  is  a  dimensionless  indicator  of  the  degree 
to  which  the  wind  waves  have  reached  steady  state  and  is 
frequently  related  to  the  non-dimensional  fetch  as  (Wu, 
1985)  : 


(CQ/U*2)  =  0.05[gx/(U*2)]3/10 


(53) 


Hsu  (1974)  related  the  wave  saturation  parameter  to  the 
roughness  length,  zQ,  by  defining  the  Charnock  coefficient 


in  terms  of  the  slope  of  the  wave  of  frequency  nQ '  .  His 


formulation  became: 


(H/2,r)(C0/U*)' 


(54) 


where  H  is  the  significant  wave  height.  zQ  can  be  seen  to 
vary  according  to  the  local  wind  wave  spectral  peak. 

Numerous  wave  models  have  incorporated  these  relation¬ 
ships  in  an  attempt  to  relate  variations  in  the  surface  wave 
energy  spectrum  to  variations  in  the  magnitude  of  the  drag 
coefficient.  Davidson  (1974)  formulated  an  empirical  model 
in  which  the  drag  coefficient  is  represented  by: 


[k/[ln(z/z0)  -  xp  *  b[(C0/U*)  -  26. 3  ]  ]  ] 2 


(55) 


where  b  is  an  empirical  constant  found  to  be  0.13. 

Others  have  employed  indirect  methods  to  calculate  Cp 
using  roughness  length  parameterizations .  Byrne  (1982) 
derived  a  roughness  length  in  the  form  of  an  integrated 
slope  spectrum: 

zQ  =  Cb/J’"  S(n'  )  (n’  )2dn'  (56) 

where  Cp  is  an  adjustable  coefficient. 

Kitaigorodskii  (1973)  derived  a  roughness  length  model 
in  the  following  form: 

z„  =  Z/mCmn')e-2ks/a'V*dn’)UZ  (57) 

One  purpose  of  the  TOWARD  Experiment  was  to  evaluate 
scatterometer  return  signals  in  order  to  examine  surface 
roughness  as  a  predictor  of  windspeed  over  the  ocean.  These 
signals  can  be  correlated  with  windspeed  since  they  increase 
with  increasing  surface  roughness.  The  surface  roughness  is 
dependent  on  r,  which  from  Eqn.  (2a)  is  a  function  of  Cp  and 
the  ten  meter  height  windspeed,  U^g.  The  scatterometer, 
however,  evaluates  U-^g  based  on  signal  reflection  from  the 
surface  (Geernaert,  1983),  while  the  drag  coefficient  has 
been  treated  as  a  constant  or  as  a  simple  function  of  U-^g  . 
Thus,  the  dependence  of  the  drag  coefficient  on  the  surface 
roughness  is  not  fully  understood.  Because  wave  growth  and 
decay  affect  the  magnitude  of  Cp  (Geernaert,  1985a),  scat¬ 
terometer  predictions  of  windspeed  based  on  the  bulk 
formulae  could  be  improved  by  including  roughness  as  an 
additional  variable. 

The  TOWARD  Experiment  was  designed,  in  part,  to  look  for 
a  trend  in  Cp  as  a  function  of  time  over  the  course  of  the 
sea  breeze.  The  sea  breeze  is  characterized  by  growing  winds 
and  sea  state  along  with  corresponding  increases  in  surface 
roughness.  By  observing  the  behavior  of  this  coefficient 
with  time,  the  interaction  between  wind  and  wave  state  may 
be  better  understood. 


III.  SYNOPTIC  METEOROLOGY  DURING  THE  TOWARD  EXPERIMENT 


The  TOWARD  experiment  was  conducted  at  the  Naval  Ocean 
Systems  Center  (NOSC)  tower  which  is  located  two  kilometers 
offshore  of  Mission  Beach,  San  Diego.  Minimal  surface  layer 
climatological  data  is  available  for  the  site.  Blanc  (1981) 
summarized  the  environmental  conditions  for  San  Nicholas 
Island  (SNI),  California,  which  is  located  roughly  150  km 
northwest  of  the  NOSC  tower.  While  wind  and  wave  conditions 
are  more  vigorous  at  SNI  than  those  at  the  tower,  the  air 
and  sea-surface  temperatures  are  consistent  between  the  two 
sites.  Blanc  determined  the  most  frequent  true  wind  direc¬ 
tion  to  be  from  the  northwest  throughout  the  year  with  mean 
windspeeds  of  7.7  m/s  in  March  and  6.2  m/s  for  both  October 
and  November.  During  the  months  of  March,  October  and 
November,  when  the  TOWARD  experiment  was  conducted,  the  mean 
climatological  daily  temperatures  (°C)  are  15.8  for  March, 
20.4  for  October  and  18.4  for  November,  with  average  daily 
minimum  temperatures  of  12.4,  16.0  and  14.9,  respectively. 

Mean  climatological  sea-surface  temperatures  for  the  same 
months  are  14.4  for  March,  with  values  of  18.2  and  16.6  for 
October  and  November. 

The  mean  conditions  for  San  Nicholas  Island  differ  some¬ 
what  from  those  reported  onshore  at  Lindbergh  Field  in  San 
Diego.  The  mean  daily  maximum  temperatures  for  March 
October  and  November  are  18.9,  23.3  and  21.2,  respectively, 
while  the  average  daily  minimum  temperatures  are  10.1,  14.7 

and  10.8.  In  March  and  October,  the  prevailing  winds  are 
generally  west-northwesterly  but  shift  to  a  northeasterly 
direction  in  November.  Mean  windspeeds  are  also  smaller 
with  values  of  3.3,  2.9  and  2.5  m/sec  reported  for  March, 
October  and  November,  respectively. 


The  TOWARD  Experiment  was  designed  to  cover  a  large 
range  of  wind  conditions,  particularly  for  high  winds. 
Measurement  periods  were  therefore  chosen  to  be  from 
mid- September  through  mid-November  1984  (Phase  I)  and  March 
1985  (Phase  II). 

The  almost  daily  presence  of  the  local  sea  breeze  is  a 
factor  in  the  environmental  conditions  at  the  NOSC  tower. 
Sea  breeze  circulations  develop  when  air  over  land  is  heated 
much  more  than  the  air  over  the  ocean.  Since  the  land 
usually  becomes  warmer  than  the  sea  during  the  day,  density 
differences  over  the  two  surfaces  bring  cool  air  from  the 
sea  toward  the  warmer  land  at  low  levels.  A  return  flow  of 
warm  air  from  the  land  to  the  ocean  occurs  at  upper  levels. 
At  night,  radiative  cooling  takes  place  over  the  land  while 
the  sea  surface  maintains  a  more  uniform  temperature  due  to 
the  large  heat  capacity  of  the  ocean.  The  result  is 
nocturnal,  low  level  flow  from  the  land  toward  the  sea 
(Hess,  1979). 

Because  the  largest  horizontal  density  gradient  occurs 
at  the  coastline,  the  sea  breeze  dominates  in  the  coastal 
zone,  while  the  geostrophic  wind  (V  )  becomes  a  more  impor- 
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tant  factor  for  describing  surface  layer  energetics  farther 
offshore.  It  follows  that  waves  generated  by  the  sea  breeze 
would  tend  to  be  of  relatively  short  fetch.  Fig.  3.11  from 
31  October  provides  an  example  of  synoptic  conditions  for 
Southern  California  during  the  TOWARD  Experiment.  The  large 
scale  flow  (Vg)  can  be  seen  to  originate  from  the  northwest. 
Ship  observations,  300-400  km  offshore,  also  indicate  a 
northwesterly  flow.  These  wind  observations  from  offshore 
sites  show  little  diurnal  variation.  At  1000  local  time, 
however,  reports  from  Lindbergh  Field,  San  Diego,  show  winds 
from  a  south- southwesterly  direction,  indicating  smaller 
scale,  local- forcing  nearer  the  coast.  A  significant  diurnal 
shift  in  wind  direction  can  also  be  seen  in  these  reports. 


Thus,  the  sea  breeze  appears  to  exert  an  influence  less  than 
300  to  400  km  from  the  shore. 

Since  the  sea  breeze  is  of  a  different  scale  than  the 
synoptic,  geostrophic  wind,  computation  of  Vg  provides  an 
indication  of  the  relative  magnitudes  of  these  two  types  of 
forcing  near  the  NOSC  tower.  Calculations  of  the  geos¬ 
trophic  wind  from  synoptic  charts  were  made  for  1000  local 
(1800  GMT)  and  1600  local  (2400  GMT). 

Vg  =  (l/fpa)(dp/dn)  (58) 

where  Vg  represents  the  geostrophic  wind,  f  is  the  coriolis 
parameter,  p  is  the  atmospheric  pressure  and  n  is  a  unit 
vector  normal  to  the  isobars.  These  times  were  chosen  in 
order  to  compare  the  differences  in  sea  breeze  magnitudes 
between  morning  and  afternoon.  Typically,  the  strongest  sea 
breeze  occurs  in  mid-afternoon  when  the  land  surface  has 
heated  to  its  daily  maximum.  Since  final  surface  charts 
were  not  available  for  2100  GMT  (1300  local),  analyses  for 
2400  GMT  were  used  instead.  These  values  were  compared  with 
windspeed  observations  collected  at  the  tower  and  with  data 
from  ship  reports  in  the  Southern  California  area. 

Data  collected  at  the  tower  indicated  a  general  increase 
in  windspeed  from  1000  to  1600  local  (1800  GMT  to  2400  GMT). 
The  geostrophic  wind  velocities  based  on  Eqn.  (58)  and  the 
synoptic  pressure  maps  were  smaller  and  showed  much  less 
diurnal  variation  than  those  of  the  tower  (Vt).  Ship  obser¬ 
vations  from  300-400  km  offshore  indicated  consistently 
higher  windspeeds  for  both  time  periods,  usually  ranging 
from  5-15  kt .  Although  the  computed  V  and  ship  reported 

O 

surface  winds  (Vs)  showed  considerable  differences  in  magni¬ 
tude,  the  wind  directions  from  offshore  sites  during  the 
experiment  showed  no  statistically  significant  diurnal 
effects.  Tables  II  and  III  show  these  results. 


TABLE  II 

Computed  Vg  &  Ship  Reported  Sfc  Winds  (Phase  I) 


(V 

t-  vg 

and  Vs 

in  m/ s ) 

Date 

Time 

vt 

Dir. 

Vg 

Dir. 

vs 

Dir. 

D. (km) 

10/2 

hi 

Z 

2.0 

87 

<1 

NW 

2.6 

N 

400 

Z 

-  -  - 

— 

1.2 

NW 

7.7 

N 

400 

10/3 

18 

Z 

2.5 

297 

1.4 

NW 

— 

— 

— 

24 

Z 

5.3 

309 

<1 

NW 

7.7 

N 

400 

10/6 

18 

Z 

-  -  - 

— 

<1 

NNW 

2.6 

W 

400 

24 

Z 

5.9 

313 

1.4 

NNW 

7.7 

N 

400 

10/7 

18 

Z 

4.6 

335 

<1 

NW 

2.6 

N 

400 

24 

Z 

— 

— 

1.7 

NW 

5.2 

N 

400 

10/11 

18 

Z 

— 

— 

1.4 

NW 

7.7 

N 

400 

24 

Z 

— 

— 

1.1 

NW 

7.7 

NW 

400 

10/12 

18 

Z 

1.0 

154 

1.0 

NE 

10.2 

N 

400 

24 

Z 

8.3 

317 

1.4 

NE 

10.2 

N 

400 

10/15 

18 

Z 

6.0 

180 

1.4 

NE 

10.2 

N 

400 

24 

Z 

5.4 

293 

1.2 

NE 

12.8 

N 

400 

10/16 

18 

Z 

— 

— 

<1 

NNW 

2.6 

NW 

400 

24 

Z 

5.4 

203 

<1 

NW 

7.7 

N 

400 

10/17 

18 

Z 

6.7 

280 

1.4 

N 

10.2 

N 

400 

24 

Z 

6.2 

181 

1.7 

N 

7.7 

NW 

350 

10/18 

18 

Z 

— 

1.0 

ENE 

10.2 

7.7 

NNW 

400 

24 

Z 

6.8 

337 

<1 

ENE 

N 

400 

10/20 

18 

Z 

5.6 

255 

<1 

N 

— 

— 

— 

24 

Z 

6.0 

267 

1.4 

NW 

10.2 

N 

400 

10/21 

18 

24 

Z 

Z 

1.8 

127 

<1 

NNE 

NW 

7~7 

N 

400 

10/23 

18 

Z 

— 

— 

<1 

NE 

5.1 

N 

400 

24 

Z 

8.7 

288 

<1 

ENE 

— 

— 

— 

10/25 

18 

z 

— 

— 

<1 

SE 

2.6 

N 

400 

24 

z 

— 

— 

<1 

SE 

2.6 

NW 

400 

10/27 

18 

z 

— 

— 

<1 

N 

— 

— 

— 

24 

z 

4.4 

221 

1.7 

N 

5.2 

N 

400 

10/28 

18 

z 

— 

— 

<1 

NE 

— 

— 

— 

24 

z 

3.9 

298 

<1 

NW 

7.7 

N 

400 

10/29 

18 

z 

— 

— 

1.1 

NE 

— 

— 

— 

24 

z 

4.4 

315 

<1 

NE 

5.2 

NW 

400 

10/30 

18 

z 

1.4 

231 

<1 

NNW 

— 

— 

— 

24 

z 

3.7 

294 

2.8 

NW 

7.7 

N 

400 

10/31 

18 

z 

2.8 

218 

1.7 

NW 

— 

— 

— 

24 

z 

5.3 

296 

1.7 

NW 

— 

— 

— 

11/2 

18 

z 

2.6 

298 

1.7 

NNW 

7 . 7 

N 

400 

24 

z 

6.8 

325 

1.4 

NW 

7.7 

N 

400 

11/3 

18 

z 

— 

— 

<1 

NW 

2.7 

SW 

400 

24 

z 

3.9 

289 

1.4 

NE 

10.2 

N 

400 

11/4 

18 

z 

1.1 

170 

<1 

SE 

7  .  7 

N 

400 

24 

z 

2.8 

306 

<1 

SE 

— 

— 

— 

11/5 

18 

z 

1.4 

158 

<1 

N 

5.1 

N 

400 

24 

z 

1.5 

199 

2.1 

NW 

7 . 7 

N 

400 

11/7 

18 

z 

— 

— 

<1 

NW 

7.7 

N 

400 

24 

z 

4.1 

257 

1.4 

NW 

5.1 

N 

400 
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TABLE  III 

Computed  V_  &  Ship  Reported  Sfc  Winds  (Phase  II) 


(V 

o  vg 

and  V„ 

in  m 

/s) 

Date 

Time 

vt 

Dir . 

Vg 

Dir. 

vs 

Dir. 

D. (km) 

3/12 

18 

Z 

•  -  • 

-  -  - 

1.4 

NNW 

5.1 

N 

400 

24 

Z 

6.0 

260 

<1 

NW 

10.8 

N 

400 

3/13 

18 

Z 

— 

- - 

<1 

NE 

7.7 

NW 

250 

5.1 

NW 

275 

24 

Z 

3.0 

260 

1.9 

N 

7.7 

N 

400 

3/15 

18 

Z 

— 

— 

<1 

NW 

5.1 

N 

400 

24 

Z 

4.8 

415 

<1 

NW 

7.7 

N 

400 

3/16 

18 

Z 

4.0 

003 

<1 

NW 

2.6 

NW 

400 

24 

Z 

4.9 

340 

<1 

NW 

7.7 

N 

400 

3/18 

18 

Z 

— 

— 

1.4 

NW 

7.7 

N 

400 

24 

Z 

9.4 

296 

2.4 

NW 

10.3 

N 

400 

3/19 

18 

Z 

— 

— 

2.8 

NW 

10.3 

N 

400 

24 

Z 

4.8 

275 

1.7 

NW 

12.9 

N 

400 

3/20 

18 

Z 

- - 

— 

1.3 

NNW 

7.7 

NW 

400 

24 

Z 

5.9 

312 

1.4 

NW 

7 . 7 

N 

400 

3/21 

18 

Z 

— 

— 

2.1 

NW 

10.3 

N 

400 

24 

Z 

3.2 

315 

1.8 

NW 

18.0 

N 

400 

3/22 

18 

Z 

— 

— 

1.4 

N 

12.9 

N 

400 

24 

Z 

3.6 

322 

1 

NE 

12.9 

N 

400 

3/24 

18 

Z 

— 

— 

<1 

NE 

5.1 

NW 

275 

5.1 

N 

400 

24 

Z 

— 

— 

<1 

NNW 

7 . 7 

N 

400 

3/25 

18 

Z 

— 

— 

1.4 

NNW 

7.7 

N 

400 

24 

Z 

6.1 

220 

1.7 

NNW 

10.3 

N 

400 

3/26 

18 

Z 

3.6 

218 

1.7 

NW 

7.7 

N 

400 

24 

z 

3.6 

265 

<1 

NW 

10.3 

N 

400 

3/27 

18 

z 

7.6 

233 

<1 

NW 

5.1 

NNW 

400 

24 

z 

5.5 

264 

2.1 

NW 

7.7 

NNW 

400 

3/28 

18 

z 

7.5 

235 

2.8 

NW 

12.9 

N 

400 

24 

z 

— 

— 

2.8 

NW 

10.3 

N 

400 

A. .  THE  TOWARD  SYNOPTIC  AND  SURFACE  LAYER  CONDITIONS 


This  section  will  be  a  presentation  of  the  surface  layer 
and  larger  scale  meteorological  conditions  for  the  TOWARD 
experiment.  The  data  set  is  divided  into  two  sections,  Phase 
I  and  Phase  II.  Phase  I  covers  the  period  from  2  October  to 
8  November  1984,  while  Phase  II  covers  the  period  of  12  to 
28  March  1985.  Each  phase  is  discussed  separately. 
Individual  surface  analyses  for  1800  GMT  will  be  presented 


in  the  text  for  days  on  which  data  were  collected. 
Additional  figures  representing  wind  speed  and  direction,  as 
well  as  air  and  sea- surface  temperatures  can  be  found  in 
Appendix  A.  The  latter  figures  present  time  series  of 
discussed  quantities  and  provide  a  reference  for  the 
following  discussion.  Unless  otherwise  noted,  the  times  are 
local  (PDT  for  Phase  I  and  PST  for  Phase  II). 

1.  Phase  I 

During  the  first  week  of  the  TOWARD  experiment, 
winds  were  from  the  northwest  and  windspeeds  were  light  to 
moderate.  A  maximum  value  of  roughly  7  m/s  (Fig.  A.l)  was 
recorded  on  October  6th.  Strong  local  diurnal  variations  due 
to  the  sea  breeze  also  prevailed  in  the  first  week,  which  is 
particularly  evident  in  the  data  from  3  and  6  October.  On 
those  dates,  windspeeds  increased  Significantly  from  0900 
PDT,  remained  steady  for  several  hours,  and  then  decreased 
markedly  throughout  the  sampling  period  which  ended  at  2215. 

Data  collection  resumed  at  1235  PDT  on  October  6th. 
Recorded  windspeeds  remained  fairly  steady  from  this  time  to 
1500,  when  values  began  decreasing.  Data  were  only  collected 
for  one  hour  on  7  October.  Windspeeds  increased  steadily 
from  about  4  to  6  m/s  over  this  period. 

Synoptic  charts  for  the  2-8  October  period  do  not 
indicate  any  significant  weather  activity.  A  low  pressure 
center  was  present  southeast  of  the  San  Diego  area  early  in 
the  week,  before  moving  farther  southeast.  Weather  condi¬ 
tions  were  subsequently  dominated  by  a  strong  offshore  ridge 
which  persisted  throughout  the  week. 

Data  were  collected  only  on  the  11th  and  12th  of 
October  during  the  second  week  of  TOWARD.  Wind  directions 
were  from  the  northwest  on  the  11th  and  windspeeds  were  low. 
On  the  12th,  winds  were  also  out  of  the  northwest,  but  had 
much  higher  velocities.  A  maximum  of  about  9  m/s  was 


Figure  3.1  Surface  analyses  1800 
and  3  (b)  October  1984. 


Figure  3.2  Surface  analyses  1800  GMT  6  (a) 
and  7  (b)  October  1984. 
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recorded  on  that  date  at  1550.  Typical  sea  breeze  patterns 
were  also  apparent  in  the  data,  particularly  on  the  12th. 
Weak  easterly  winds  were  present  in  the  early  morning.  The 
winds  shifted  to  the  northwest  near  noon  and  increased 
rapidly  to  a  maximum.  These  values  began  to  drop  off  at 
approximately  1600  PDT. 

The  synoptic  charts  from  this  period  show  a  frontal 
passage  through  the  San  Diego  area  on  the  morning  of  12 
October.  This  event  cannot  be  correlated  with  the  strong 
winds  of  the  12th,  because  the  front  was  already  several 
hundred  kilometers  to  the  southeast  when  the  winds  were 
recorded.  Synoptic  charts  do  not  indicate  the  presence  of  a 
strong,  large-scale  pressure  gradient  associated  with  the 
trough.  Thus  the  most  likely  cause  of  these  high  winds  was 
strong  diurnal  heating  and  subsidence  following  the  passage 
of  the  front. 

Data  collection  for  the  third  week  of  the  experiment 
began  at  0900  PDT  on  October  15th.  A  definite  wind  shift  can 
be  seen  on  the  15th  accompanied  by  sharp  fluctuations  in 
windspeed.  Synoptic  charts  for  that  date  indicated  the  pres¬ 
ence  of  a  stationary  front  almost  directly  over  the  San 
Diego  area  at  1000  PDT.  This  front  passed  through  the  area 
at  a  rate  of  only  2  m/s.  Its  slow  rate  of  movement  caused  it 
to  influence  weather  conditions  in  the  San  Diego  area 
through  the  16th. 

Wind  shifts  plus  strong  increases  and  decreases  in 
windspeed  were  also  apparent  on  the  16th  and  17th.  Synoptic 
charts  for  these  dates  show  a  front  moving  rapidly  through 
the  tower  area  at  a  rate  of  over  10  m/s.  This  front  was  well 
to  the  southeast  within  six  hours. 

October  18th  provided  a  more  typical  weather  pattern 
with  northwesterly  winds  and  a  strong  increase  in  windspeed 
between  1000  and  1700  local  due  to  the  diurnal  heating 
pattern.  Synoptic  charts  show  a  strong  high  pressure  ridge 


Figure  3.3  Surface  analyses  1800  GMT  11  (a) 
and  12  (b)  October  1984. 


building  to  the  southwest,  with  a  weaker  ridge  to  the 
northeast . 

Winds  were  westerly  on  the  morning  of  the  20th  with 
a  shift  to  an  easterly  flow  near  midnight.  Velocities  were 
moderate  and  fairly  constant,  with  values  decreasing  in  late 
afternoon. 

The  first  observations  for  the  week  of  22-29  October 
began  at  1100  on  October  23rd.  A  strong  increase  in  wind- 
speed  can  be  seen,  along  with  a  northwesterly  wind  direc¬ 
tion.  On  the  afternoon  of  the  25th,  windspeeds  again  showed 
an  increase.  The  wind  direction  was  southerly  in  the  late 
morning  and  shifted  to  the  northwest  by  early  afternoon. 
Wind  speeds  recorded  on  the  27th  were  moderate  with  little 
diurnal  change.  This  pattern  continued  on  the  28th,  although 
the  winds  again  originated  from  the  northwest. 

The  fifth  week  of  the  TOWARD  experiment  was  from  29 
October  to  5  November.  Windspeeds  were  moderate  on  the  29th 
and  30th  and  were  from  the  northwest.  A  wind  shift  and 
increase  in  windspeed  began  at  approximately  0900  PDT  on  the 
31st.  Evaluation  of  synoptic  charts  for  31  October  indicated 
the  presence  of  a  large-scale  pressure  gradient  associated 
with  an  oncoming  front.  On  November  2nd,  the  winds  were 
again  out  of  the  northwest  and  San  Diego  returned  to  the 
typical  sea  breeze  pattern. 

The  development  of  a  strong  high  pressure  center 
over  the  western  United  States  produced  the  low  windspeeds 
and  shift  in  direction  evident  on  4  and  5  November. 
Windspeeds  ranged  from  approximately  2  to  6  m/s. 

Data  collection  for  the  final  week  of  Phase  I 
concluded  on  7  November.  Winds  on  the  5th  were  light  and 
southwesterly.  By  the  7th,  the  typical  sea  breeze  pattern 
again  prevailed.  Windspeeds  were  low  and  out  of  the  west  in 
early  morning,  but  shifted  to  the  northwest  by  early  after¬ 
noon.  A  maximum  value  of  about  5  m/s  was  recorded  at  1105 
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Figure  3.7  Surface  analyses  1800  GMT  23  (a) 
and  25  (b)  October  1984. 


Figure  3.10  Surface  analyses  1800  GMT  31  (a) 
and  2  (b)  November  1984. 
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Figure  3.11  Surface  analyses  1800  GMT  3  (a) 
and  4  (b)  November  1984. 


48 


2.  Phase  II 


Data  collection  for  Phase  II  of  the  TOWARD  experi¬ 
ment  began  on  12  March  1985.  It  can  be  seen  from  Fig.  3.14 
that  winds  on  the  12th  and  13th  were  light  to  moderate  and 
were  primarily  westerly.  On  March  15th  and  16th  winds  were 
again  moderate,  but  had  a  more  typical  northwesterly  flow. 

Synoptic  charts  for  this  period  do  not  indicate  any 
significant  weather  activity.  Conditions  were  dominated  by 
an  offshore  high  pressure  center  throughout  the  week. 

On  March  18th,  synoptic  charts  depict  a  tight  pres¬ 
sure  gradient  resulting  from  the  passage  of  a  low  pressure 
center  through  the  San  Diego  area.  The  result  was  a  sharp 
increase  in  windspeed.  A  maximum  value  of  roughly  10  m/s  was 
recorded  at  1550  on  that  date.  During  the  remainder  of  the 
week,  the  development  of  a  strong  offshore  ridge  reduced 
windspeeds  and  produced  a  northwesterly  flow. 

The  final  week  of  the  TOWARD  experiment  covered  25 
to  28  March.  A  gradual  wind  shift  accompanied  by  light  to 
moderate  winds  can  be  seen  on  the  25th.  Synoptic  charts  for 
that  date  show  a  front  passing  through  the  San  Diego  area 
between  1000  and  1600  local  time,  at  an  approximate  rate  of 
4-5  m/ s .  As  it  moved  westward,  the  low  deepened  and  the 
frontal  band  remained  southeast  of  San  Diego  through  1000  on 
the  26th.  The  winds  reflect  this  frontal  movement.  A  shift 
in  direction  to  a  more  typical  northwesterly  flow  occurs  in 
late  morning  and  a  sea  breeze  pattern  is  apparent  by  late 
afternoon.  The  maximum  wind  speed  recorded  was  about  7  m/s 
at  1715. 

Winds  on  the  27th  were  out  of  the  west  -  southwest  and 
had  generally  higher  magnitudes  that  those  recorded  the 
previous  day.  These  values  ranged  from  a  mimimum  of  about  3 
m/s  to  a  maximum  of  just  over  8  m/s  at  1015.  This  maximum  is 
supported  by  synoptic  charts  for  1000  that  depict  a  strong, 
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large-scale  pressure  gradient  associated  with  a  deep  low 
center  over  Vancouver,  Canada.  This  system  quickly  moved 
southeast,  where  the  low  center  deepened  over  Arizona  and 
New  Mexico  for  the  next  several  days.  The  result  was  a  tight 
pressure  gradient  which  produced  the  high  wind  values 
collected  on  28  March.  On  that  date,  a  maximum  wind  value  of 
over  12  m/ s  was  recorded  at  1415. 


Figure  3.18  Surface  analyses  1800  GMT  25  (a) 
and  26  (b)  March  1985. 


IV.  DRAG  COEFFICIENT  RESULTS 


A  principal  motivation  of  this  thesis  is  to  examine  the 
variation  of  the  sea-surface  drag  coefficient  with  varying 
wind  and  wave  states.  Since  windspeed  magnitudes  were  gener¬ 
ally  low  to  moderate  (i.e.  <  10-11  m/s)  during  both  phases 
of  TOWARD,  the  primary  emphasis  will  be  placed  on  CgN  varia¬ 
tions  as  a  function  of  wave  state. 

The  question  of  whether  the  neutral  drag  coefficient  is 
constant  or  increases  with  windspeed,  wave  state  and  fetch 
has  received  considerable  study  in  recent  years.  Garratt 
(1977)  summarized  sea-surface  drag  coefficients  from  17 
experiments  (based  on  eddy  correlation  and  profile  methods) 
published  between  1967  and  1975.  Observations  of  Cj-^  (repre¬ 
sentative  of  a  height  of  10  m)  plotted  as  a  function  of  U-^g 
showed  that  Charnock’s  relation  (Eqn.  26)  was  valid  over  a 
windspeed  range  of  4  to  21  m/s. 

Large  and  Pond  (1981)  compared  measurements  of  the 
momentum  flux  obtained  by  both  Reynolds  flux  and  dissipation 
methods.  They  reported  that  a  comparison  of  both  methods 
showed  excellent  agreement,  on  average,  for  windspeeds  of 
4-20  m/ s .  Large  and  Pond  also  plotted  neutral  drag  coeffi¬ 
cients  against  U^q.  Although  they  found  a  slight  increase  in 
the  average  CDN  with  windspeeds  below  10  m/s,  they  felt  that 
this  relationship  could  be  adequately  described  by  a 
constant,  where: 

103Cdn  =  1.14  4<U10<10  m/s  (59) 

Above  10  m/s,  the  values  of  CDN  demonstrated  a  more  rapid, 
linear  increase  up  to  windspeeds  of  26  m/s. 

Smith  (1980)  used  the  eddy  correlation  method  to  obtain 
windspeed  and  heat  flux  measurements  over  the  windspeed 
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range  of  6-22  m/s.  He  noted  that  the  Chamock  relation  also 
held  for  neutral  sea-surface  drag  coefficients  measured  for 
long  fetch.  These  Cpjj  values  were  clearly  seen  to  increase 
with  increasing  windspeed.  Several  points  were  not 
adequately  explained  by  the  Charnock  hypothesis,  however. 
These  included: 

1.  The  dependence  of  the  drag  coefficient  on  windspeed 
was  greater  than  predicted  if  the  data  at  low  windspeeds 
were  not  included. 

2.  If  the  data  in  Garratt's  summary  are  viewed  sepa¬ 
rately  at  windspeeds  both  above  and  below  10  m/s,  the  drag 
coefficient  also  appears  to  be  nearly  constant  at  windspeeds 
below  10  m/ s .  Smith's  average  CDN  value  for  low  windspeeds 

*3 

was  (1.11)10  .  At  higher  windspeeds,  this  coeffcient  tends 
to  increase  rapidly  as  a  function  of  windspeed. 

Thus,  all  three  studies  indicate  a  relatively  constant 
Cjjn  value  at  windspeeds  below  10  m/s.  It  follows  then  that 
variations  in  the  neutral  drag  coefficient  at  low  to 
moderate  windspeeds  must  be  due  to  changes  in  the  sea  state. 

Garratt  did  not  explore  the  behavior  of  CDN  with  respect 
to  fetch  or  sea  state.  He  felt  that  current  quantitative 
information  was  insufficient  to  adequately  determine  the 
dependence  of  the  drag  coefficient  on  factors  other  than 
windspeed . 

Large  and  Pond  found  the  to  be  largely  independent 
of  fetch.  For  both  long  and  short  fetches,  the  mean  drag 
coefficient  remained  at  1.14  x  10  for  windspeeds  less  than 
10  m/s. 

Smith's  comparison  of  with  fetch  showed  higher  drag 
coefficients  at  short  fetches  than  those  at  longer  fetches. 
This  is  to  be  expected  since  at  short  fetches  actively 
growing  waves  absorb  momentum  from  the  wind  while  waves  at 
longer  fetches  tend  to  be  more  in  equilibrium  with  the  wind 


Smith  also  considered  the  variation  of  the  neutral  drag 
coefficient  with  sea  state.  A  linear  regression  of  Cj^q  vs. 
rms-wave  height  showed  an  increase  in  CDjj  with  an  increase 
in  wave  height,  although  the  correlation  coefficient  was 
rather  low  ( .  58 ) . 

More  recent  studies  have  included  a  drag  coefficient 
model  discussed  by  Geemaert  (1985a).  Geemaert  found  the 
model  exhibited  considerable  sensitivity  to  fetch.  The 
influence  of  fetch  on  the  model  calculations  was  due  to  the 
wave  spectral  variations  incurred  by  the  dependence  of  the 
Phillips  coefficient  on  the  wave  saturation  parameter, 
Co/U*,  where  Co  represents  the  phase  speed.  Geernaert 
observed  that  high  drag  coefficients  often  resulted  from 
fetch  limited  conditions. 

Wu  (1985)  also  looked  at  fetch  as  influenced  by  the 
Phillips  spectral  coefficient.  He  found  larger  spectral 
coefficients  at  shorter  fetches  than  those  at  longer 
fetches . 

Measurements  of  the  neutral  drag  coefficient  over  the 
sea  by  all  methods  often  yield  considerable  scatter  and  the 
data  must  be  carefully  edited.  Several  types  of  data  have 
been  removed  from  this  analysis  based  on  dissipation 
methods.  They  include: 

1.  Data  collected  during  rain  or  fog  conditions.  Water 
droplets  striking  the  hot  film  result  in  latent  heat 
release.  This  sudden  cooling  induces  a  large  jump  ^n 
voltage,  producing  spikes  in  the  data.  In  the  inertial 
subrange  particularly,  noise  from  droplet  contact  contami¬ 
nates  the  high  frequencies  and  results  in  very  large  drag 
coefficient  estimates. 

2.  Periods  of  light  winds  also  produce  erratic 
values  as  a  result  of  convection  off  the  tower  mast. 

3.  Flow  distortion  is  a  problem  with  light  or  moderate 
winds,  particularly  from  the  south.  The  best  direction  for 
data  collection  is  from  west  to  north. 


TABLE  IV 

CDN  vs.  U]_g  for  Long  and  Short  Fetches 
SHORT  FETCH  CONDITIONS :  360-115  DEGREES 


1. 0-2.0  m/s 

2.0-3 .0  m/s 

3 . 0-4 . 0  m/ s 

(1  point) 

(8  points) 

(21  points) 

value=5 . 027 

mean=2 . 017 

mean= 1.656 

s . d. = - 

s . d . = . 85 1 

s . d . = . 407 

LONG 

FETCH  CONDITIONS:  250-330 

DEGREES 

1. 0-2.0  m/s 

2 . 0- 3 . 0  m/ s 

3 . 0-4 . 0  m/ s 

(1  point) 

(53  points) 

(81  points) 

value=2 . 119 

mean: 1.485 

mean= 1.432 

s  .  d .  = - 

s . d. = . 828 

s . d . = . 850 

4 . 0-5 . 0  m/s 

5 . 0-6 . 0  m/s 

6 . 0- 7 . 0  m/ s 

(79  points) 

(63  points) 

(26  points) 

mean: .955 

mean= . 980 

mean=.896  i 

s .d. = .442 

s . d . = . 444 

s . d .  =  . 3  77 

7 . 0- 8 . 0  m/ s 

8 . 0- 9 . 0  m/ s 

9.0-10.0  m/ s  i 

(14  points) 

(8  points) 

(4  points) 

mean= .856 

mean= .939 

mean: 1.413 

s.d.= .463 

s . d . = . 254 

10.0-11.0  m/s 

(1  point) 

value= 1.479 
s  .  d .  = - 

s . d . = . 104 

I 

1 

1 

t 

1 

fairly  high  for  the  1-4  m/s  bands,  dropped  off  between  4-8 
m/s,  and  subsequently  increased  between  8-11  m/s. 
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Overall,  the  dependence  of  the  drag  coefficient  on  wind- 
speed  was  found  to  be  insignificant  above  4  m/s,  based  on 
scatter  in  the  data  as  well  as  the  large  standard  deviation 
at  each  windspeed  interval.  Below  4  m/s,  each  there  was  an 
apparent  increase  in  the  CDN  with  a  decrease  in  windspeed. 
These  values  are  somewhat  suspect,  however,  due  to  the 
uncertainty  in  the  quality  of  the  data  at  very  low 
velocities . 

Daily  time  series  of  wind  stress  were  also  examined, 
with  both  CDN  and  U*  plotted  against  time.  Lack  of  data,  low 
winds  or  strong  stability  prevented  a  close  examination  of 
the  wind  stress  in  some  instances.  Several  case  studies  are 
described  below,  however.  Appendix  A  contains  figures  rele¬ 
vant  to  the  discussion. 

1.  Frontal  Cases:  15,  16-18  October  and  26-29  March 

Fig.  A. 6  shows  the  results  when  a  front  was  passing 
the  tower  on  15  October.  CDN  values  ranged  from  .584  to 
2.058.  Magnitudes  can  be  seen  to  rise  with  rising  winds 
until  0900,  fall  with  decreasing  winds  between  0900  to  1030, 
and  remain  fairly  steady  with  the  winds  near  constant 
values.  A  sharp  drop  in  windspeed  at  1200  did  not  produce  a 
corresponding  drop  in  Cpjj  values.  increased  during  this 
period  of  low  winds  to  a  maximum  at  1500,  and  then  dropped 
off  sharply.  Synoptically ,  Fig.  A.  6  shows  a  wind  shift 
between  1400-1500  for  southerly  to  westerly  winds  with  a 
corresponding  increase  in  windspeed.  values  began  rising 
shortly  afterward  and  continued  to  increase  up  to  the  end  of 
the  sampling  period  depite  a  drop  in  windspeed  beginning  at 
approximately  1700. 

A  similar  frontal  case  is  presented  in  Fig.  A. 7  for 
16-18  October.  Both  U*  and  can  be  clearly  seen  to 
increase  with  rising  winds  to  a  maximum  at  0600  on  the  18th. 
A  wind  shift  from  southerly  to  northwesterly  winds  occurred 


shortly  after  6  am  and  several  sharp  drops  in  windspeed  took 
place  before  the  wind  leveled  off  at  approximately  6  m/s. 
The  drag  coefficient  shows  a  similar  trend  but  does  not 
reflect  the  rapid  decreases  evident  in  the  windspeed  between 
0500  to  0800. 

On  26-29  March,  passage  of  a  front  as  well  as  south¬ 
erly  winds  produce  high  Cp^  values  with  considerable  scatter 
(Fig.  A. 31). 

In  each  of  these  frontal  cases,  Cpjj  increased  just 
before  passage  of  the  front.  This  trend  is  consistent  with 
an  earlier  study  by  Denman  and  Miyake  (1973)  in  which  drag 
coefficients  were  seen  to  increase  on  the  leading  edge  of  a 
storm,  and  then  decrease  slightly  or  remain  constant.  They 
attributed  this  behavior  to  the  nature  of  the  wave  field. 
Geernaert  (1985b)  also  observed  higher  drag  coefficients 
preceding  an  approaching  front.  He  postulated  that  wave 
fields  far  ahead  of  and  behind  cold  fronts  are  semi¬ 
independent  of  one  another  and  are  generally  in  equilibrim 
with  local  winds.  Waves  generated  behind  the  front,  often  in 
a  direction  perpendicular  to  the  baroclinic  zone,  can  under 
certain  circumstances  travel  through  the  front  and  interact 
with  the  wave  field  just  ahead  of  it.  Since  the  waves  in  the 
warm  sector  often  run  parallel  to  the  local  winds,  the 
convergence  of  the  two  fields  produces  a  much  more  energetic 
sea  state  just  ahead  of  the  front.  The  result  would  be 
higher  Cp^  estimates. 

2.  Sea  Breeze  Cases:  18,  27,  29  October  and  13,  15,  18 

March 

The  sea  breeze  effect  is  apparent  in  much  of  the 
data  from  the  TOWARD  experiment.  On  each  of  the  above  dates 
(Figs.  A. 8,  A. 12,  A. 14,  A. 22,  A. 23  and  A. 25),  the  Cp^  rises 
sharply  and  then  decreases  slowly  to  near  equilibrium. 
Described  in  terms  of  surface  roughness,  drag  coefficient 


values  are  known  to  increase  with  growing  wave  state.  As  the 
wave  field  becomes  increasingly  saturated  with  momentum, 
absorption  of  energy  slows,  resulting  in  drag  coefficients 
that  approach  a  smaller  and  more  constant  value. 

3.  Shift  from  Long  to  Short  Fetch  Conditions:  3-5 

November 

Data  from  3-5  November  depict  a  shift  from  lonr  to 
short  fetch  conditions.  Long  fetch  conditions  are  present 
until  approximately  2300  on  3  November,  as  indicated  by 
northwesterly  winds.  The  drag  coefficient  estimates  (Fig. 
A.  18)  are  correspondingly  low.  With  the  shift  to  northerly 
winds  at  2300,  a  sharp  increase  can  be  seen  in  CDN .  Although 
these  values  decline  over  the  next  three  hours,  they  again 
increase  due  to  the  variable  wind  direction  produced  by  a 
series  of  wind  shifts  on  the  4th. 

4.  Variable  Winds:  7  November 

High  CDN  values  as  result  of  variable  wind  direc¬ 
tions  can  be  seen  between  0900  to  1200  on  7  November  (Fig. 
A. 20).  Magnitudes  become  more  constant  as  winds  shift  to  a 
steadier  direction. 

5.  Long  Fetch  and  a  Cqjj  Decrease  with  Time:  16,19,  22 

March 

These  examples  indicate  long  fetch  conditions  due  to 
the  northwesterly  winds,  as  well  as  a  general  decrease  in 
C'DN  with  time  (Figs.  A. 24,  A. 26  and  A. 29).  Since  long  fetch 
waves  are  known  to  be  more  steady  state,  they  tend  to 
absorb  less  momentum  than  would  short  fetch  waves. 
Consequently,  the  drag  coefficient  estimates  have  lower 
magnitudes  and  tend  to  approach  a  constant  value  more 
rapidly  than  those  seen  in  earlier  examples. 


In  summary,  case  studies  suggest  a  dependence  on 
sea  state  as  illustrated  by  both  surface  roughness  and  the 
degree  of  equilibrium  between  U^q  and  the  waves  (i.e.  the 
amount  of  momentum  entering  the  wave  field  for  growth). 
Basic  to  this  hypothesis  on  wind-wave  coupling  is  that 
is  generally  independent  of  U^q  for  windspeeds  between  4-9 
m/ s ,  and  that  the  large  amount  of  scatter  is  associated  with 
systematic  trends  in  sea  state,  such  as  frontal  and  sea 
breeze  conditions.  This  scatter  should  be  reduced  once  wave 
state  is  incorporated  as  an  additional  parameter  in  CDN 
formulations . 


V.  SUMMARY  AND  CONCLUSIONS 


Surface  layer  measurements  covering  a  wide  range  of 
windpeeds  have  been  analyzed  for  both  phases  of  the  TOWARD 
experiment.  These  measurements  have  been  correlated  with 
relevant  mesoscale  and  synoptic  features.  The  results  were 
then  interpreted  relative  to  the  C-^. 

Neutral  drag  coefficient  values  for  both  phases  of 
TOWARD  showed  no  significant  windspeed  dependence  in  the 
range  between  4  and  9  m/s.  The  mean  value  was  determined  to 
be  (.94  +,  0.4)10  ,  which  is  somewhat  smaller  than  results 
obtained  by  Smith  (1980)  and  Large  and  Pond  (1981).  These 
differences  appear  to  be  insignificant,  however,  due  to  the 
scatter  present  in  all  three  data  sets.  CDN  values  are 
higher  for  windspeeds  less  than  4  m/sec,  but  these  values 
are  suspect  as  low  windspeeds  often  produce  unreliable  drag 
coefficient  estimates. 

Definite  trends  were  observed  in  CDjj  with  fetch  and  sea 
state.  Short  fetch  conditions  produced  generally  higher  Cj-,^ 
values  than  those  for  long  fetch.  Trends  in  sea  state  such 
as  frontal  and  sea  breeze  conditions  were  reflected  in  drag 
coefficient  estimates  as  well.  As  with  Large  and  Pond,  CDN 
estimates  were  seen  to  increase  sharply  just  prior  to  the 
passage  of  a  front.  During  sea  breeze  conditions,  the 
values  increased  with  time.  As  windspeeds  and  the  wave  field 
reached  equilibrium,  however,  the  drag  coefficients 
decreased  toward  a  more  constant  value. 
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APPENDIX  A 

TOWARD  DAILY  TIME  SERIES 


This  appendix  contains  daily  time  series  for  wind  speed 
and  direction,  air  and  sea  surface  temperatures,  the  neutral 
drag  coefficient  and  U*,  for  both  phases  of  the  TOWARD 
experiment.  These  figures  provide  a  reference  for  the 
discussion  provided  in  chapters  three  and  four. 
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